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With the possibility of testing massive gravity in the context of black hole physics in mind, we 
consider the radiation produced by a particle plunging from slightly below the innermost stable 
circular orbit into a Schwarzschild black hole. In order to circumvent the difficulties associated 
with black hole perturbation theory in massive gravity, we use a toy model in which we replace the 
graviton field with a massive scalar field and consider a linear coupling between the particle and this 
field. We compute the waveform generated by the plunging particle and study its spectral content. 
This permits us to highlight and interpret some important effects occurring in the plunge regime 
which are not present for massless fields such as (i) the decreasing and vanishing, as the mass 
parameter increases, of the signal amplitude generated when the particle moves on quasicircular 
orbits near the innermost stable circular orbit and (ii) in addition to the excitation of the quasinormal 
modes, the excitation of the quasibound states of the black hole. 

PACS numbers: 04.30.Db, 04.25.Nx, 04.70.Bw 


I. INTRODUCTION 

In recent years, theories of massive gravity, i.e. ex¬ 
tensions of general relativity mediated by an ultralight 
massive graviton, are the subject of intense activity and 
the internal inconsistencies present in the old Fierz-Pauli 
theory [1, 2] (in particular, the discontinuity with general 
relativity in the limit where the mass of the graviton is 
taken to zero and the ghost-problem) seem to be over¬ 
come (see, e.g., Refs. [3] and [4] for reviews on massive 
gravity and Ref. [5] for experimental constraints on the 
mass of the graviton). Because massive gravity could ex¬ 
plain naturally the observed accelerated expansion of the 
present Universe, i.e. without requiring dark energy or 
introducing a cosmological constant (see, e.g.. Sec. 12 of 
Ref. [4] and references therein), it has led to numerous 
works in cosmology. By contrast, it has rarely been the 
subject of works in black hole (BH) physics (see, however. 
Ref. [6] and references therein for recent articles dealing 
with BH solutions in massive gravity and Refs. [7-9] for 
important considerations on the problem of BH stability 
in massive gravity). In our opinion, it is in this par¬ 
ticular area that, in the near future, theories of massive 
gravity could be directly tested with the next generations 
of gravitational wave detectors. 

In this article, with the possibility of testing massive 
gravity in the context of BH physics in mind, we intend to 
consider the radiation produced by a “particle” plunging 
from slightly below the innermost stable circular orbit 
(ISCO) into a Schwarzschild BH. We shall assume an ex¬ 
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treme mass ratio for this physical system, i.e. that the 
BH is much heavier than the particle. Such a hypoth¬ 
esis permits us to describe the emitted radiation in the 
framework of BH perturbations. Here, two remarks are 
in order: 

— In the context of Einstein’s general relativity, the 
problem we consider is of fundamental importance 
and there exists a large literature concerning it 
more or less directly (see, e.g.. Refs. [10-25]). In¬ 
deed, the “plunge regime” is the last phase of the 
evolution of a stellar mass object orbiting near a 
supermassive BH (note that it can also be used to 
describe the late-time evolution of a binary BH) 
and the waveform generated during this regime en¬ 
codes the final BH fingerprint. 

— The Schwarzschild BH, one of the most important 
solutions of Einstein’s general relativity, is also rel¬ 
evant to massive gravity [6]. Indeed, it is a solu¬ 
tion of the pathology-free bimetric theory of grav¬ 
ity discussed by Hassan, Schmidt-May and von 
Strauss in Ref. [26] and obtained by extending, 
in curved spacetime, the fundamental work of de 
Rham, Gabadadze and Tolley [27, 28]. 

In the context of the bimetric theory of gravity dis¬ 
cussed in Ref. [26], gravitational perturbations of the 
Schwarzschild BH have been studied by Brito, Cardoso 
and Pani [8, 9]. This problem is rather difficult and leads 
to results which generalize the problem of massless spin-2 
fluctuations in a nontrivial way. In particular, it is im¬ 
portant to recall that (i) the Schwarzschild BH is linearly 
unstable for small graviton masses, (ii) the € = I modes 
(the so-called dipole modes) become dynamic and (iii) 
except for the odd-parity dipole modes, all the other dy- 
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namical modes are governed by two or three coupled dif¬ 
ferential equations. All these results make very difficult, 
in the context of massive gravity, the theoretical analy¬ 
sis of the gravitational radiation produced by a particle 
falling into a Schwarzschild BH on an arbitrary trajec¬ 
tory. To our knowledge, there does not exist any work 
dealing fully with this fundamental problem. In a re¬ 
cent paper [29] (see also Ref. [30] for a more complete 
study including analytical results and extension to other 
bosonic fields), we have partially considered it by only 
focusing on some aspects linked with the excitation of 
quasinormal modes (QNMs). We have shown that, in 
rather large domains around critical values of the gravi¬ 
ton mass, the excitation factors of the long-lived QNMs 
have a strong resonant behavior which induces, as a con¬ 
sequence, the existence of giant and slowly decaying ring¬ 
ings. This unexpected effect was obtained in a rather 
restricted context: (i) we only focused on the odd-parity 
dipole modes; (ii) the source of the distortion of the BH 
was described by an initial value problem with Gaussian 
initial data; and (iii) we did not take into account the 
full response of the BH but only the part associated with 
QNMs and therefore we neglected various important ef¬ 
fects linked with quasibound states (QBSs) which could 
blur the QNM contribution. 

Preliminary unpublished investigations have permitted 
us to understand that, in massive gravity, in addition 
to the theoretical difficulties previously mentioned, the 
problem of the excitation of a Schwarzschild BH by a 
plunging particle is plagued by numerical instabilities. 
So, in order to simplify our task, both numerically and 
theoretically, we shall consider in this article a toy model 
in which we replace the massive spin -2 perturbations with 
a massive scalar field and consider a linear coupling be¬ 
tween the particle and this field. We shall compute the 
waveform generated by the plunging particle and study 
its spectral content. This shall permit us to describe the 
excitation of the QNMs as well as that of the QBSs of the 
BH but also to show the influence of the mass on the am¬ 
plitude of the emitted signal and, more particularly, on 
the part generated when the particle moves on quasicircu¬ 
lar orbits near the ISCO. Of course, we hope to return to 
the same problem in the context of massive gravity but 
we think that the phenomena highlighted with the toy 
model are still present in this more physical and interest¬ 
ing problem. Furthermore, the lessons we have learned 
by working with the toy model should allow us to avoid 
some difficulties due to numerical instabilities. 

Our paper is organized as follows. In Sec. H, we de¬ 
scribe our toy model, introduce our notations and estab¬ 
lish theoretically the expression of the waveform emitted 
by a scalar point particle moving on an arbitrary geodesic 
by using the standard Green’s function techniques. These 
results will be useful both in Sec. HI and Appendix A. 
In Sec. HI, from the closed-form expression of the plunge 
trajectory, we first obtain the general expression of the 
waveform generated by the plunging particle and then we 
focus more particularly on the {£ = 2, m = 2) mode of 


the scalar field. We study numerically the corresponding 
waveform (the so-called quadrupolar waveform) and its 
behavior as the mass parameter increases and as the ob¬ 
server location changes and we highlight the excitation 
of the QNMs and QBSs of the BH. In the Gonclusion, we 
summarize the main results obtained in this article, we 
briefly deal with arbitrary {£, m) modes of the scalar field 
and consider possible extensions of our work. In three ap¬ 
pendixes, we gather some important results which per¬ 
mit us to analyze and interpret the numerical results of 
Sec. HI. In Appendix A, we provide a simple closed form 
for the emitted waveform when the particle lives on the 
ISCO. In Appendix B, we focus on the weakest damped 
QNM with £ = 2 and construct numerically its response 
when it is excited by the plunging particle. In Appendix 
C, we consider the first QBSs with £ = 2 and determine 
their complex frequencies. 

Throughout this article, we adopt units such that h = 
c = G = 1 and we introduce the reduced mass parame¬ 
ter (a dimensionless coupling constant) a = 2Mfi/mp'^ 
where M, /r and rup denote, respectively, the mass of the 
BH, the mass of the scalar field and the Planck mass. 


II. THE GENERAL SOLUTION OF THE 
SCALAR WAVE EQUATION WITH SOURCE 


A. Our model 


We shall consider a point particle with mass mg and 
scalar charge q moving along a world line 7 (A) in the 
Schwarzschild spacetime {M^gap)- A is an affine param¬ 
eter and {M, gap) is the exterior of the Schwarzschild BH 
of mass M defined by the metric 
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dr^ +r^da2 ( 1 ) 


where da^ denotes the metric on the unit 2 -sphere and 
with the Schwarzschild coordinates (t, r) which satisfy 
t G] — oo,-|-oo[ and r g] 2M,-|-oo[. In the following, we 
shall also use the so-called tortoise coordinate r* g] — 
00 , -|-oo[ defined from the radial Schwarzschild coordinate 
r by dr/dr^ = (I — 2M/r) and given by r*(r) = r + 
2Mln[r/(2M) —1]. We recall that the function r* =r*(r) 
provides a bijection from ] 2 M, -|-oo[ to ] — 00 , -|-oo[. 

The particle is coupled to a scalar field $ with mass g, 
and the dynamics of the system field-particle is defined 
by the action (see also Ref. [31]) 


S — iS’field “t” ‘5'part “t” 5’int 


( 2 ) 


with 
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and 


C. Construction of the Green’s function 


'5^part — 


—Too / dr 

J y 


= -mo 


dz°‘{X) dz^{\) 

-9o.p{z{X))—^—^d\ (4) 


dX 


dX 


and 


Sint= d'^x ^-g{x) p{x)^{x). 
J M 


( 5 ) 


Here Soeid = ■S'seid [4", gajs] denotes the action of the scalar 
field $. iSpart = 5'part[-2(A), ^a/?] IS the actlon of the par¬ 
ticle, r denotes its proper time and the equations z°‘ = 
z“(A) describe its world line. S'int = 5'int[-z(A), 4>,ga/j] is 
an interaction term describing the coupling of the parti¬ 
cle and the field. It should be noted that they are coupled 
via the charge density 

S‘^{x — z{t)) 


p{x) = q f dT5'^{x, z{t)) = q f dr- 

J -y J ‘y 


V-9{x) 


( 6 ) 


associated with the charged scalar particle. 


B. Wave equation and Regge-Wheeler equation 
with source 


The wave equation governing the scalar field <i> is ob¬ 
tained by extremization of the action (2). We have 

(□-/r2)$ = -p, (7) 

where the source p is given by (6). 

Due to both the staticity and the spherical symmetry 
of the Schwarzschild background, we can solve Eq. (7) by 
introducing the decompositions 

1 /•+°° 

(t, r, 9, (p) = , _ / d(jj ( 


27r J —oo 
+ 00 -\-i 

X 

i—0 m— — 


- ytm{0,p) 


( 8 ) 


and 


p{t,r,9,p) = 
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^+oo 


xJ^'K J —oo 
+00 -\-£ 

X 

i—Q rn——£ 


du) ( 


E E 


( 9 ) 


where Yf^rn{d-,'p) are the standard spherical harmonics 
with £ G N and m = —£, —£ -I- 1,..., +(.. The wave equa¬ 
tion (7) reduces then to the Regge-Wheeler equation with 
source 




— { 1 j Pijlm- (19) 


In Eq. (10), Vf(r) denotes the Regge-Wheeler potential 
given by 


V,{r) = ( I - — 
r 


2 £(£+1) 2M 


( 11 ) 


In order to solve the Regge-Wheeler equation (10), we 
shall use the machinery of Green’s functions (see Ref. [32] 
for generalities on this topic and, e.g. Ref. [33] for its use 
in the context of BH physics). We consider the Green’s 
function Gt,j^(r*,r() defined by 




Ga;r(r*,r() =-d(r* - r() (12) 


which can be written as 

1 [C^(g)CK’’(.)> 




Wi{uj) U“^(r*)C^(r(), r*>r(.. 


(13) 


Here We{uj) denotes the Wronskian of the functions 
and These two functions are linearly independent 

solutions of the homogenous Regge-Wheeler equation 


d^ 

drl 


-f - Vi{r) 


(t>uj£ — 0 . 


(14) 


When Im(w) > 0, is uniquely defined by its purely 
ingoing behavior at the event horizon r = 2M (i.e., for 
r* —>• —oo) 


while, at spatial infinity r —>■ -|-oo (i.e., for r* 
has an asymptotic behavior of the form 


(15a) 
-boo), it 


C^(’') 

X (h 


r*-* + 
(- 


UJ 


1/2 


{u))e~dPGy, + W^X^ ln(r/M)] 


+H 


(+)(a;)e 


+’i[p(a;)r* + [M/x^/p(u;)] ln(r/M)] 


) •(15b) 


Similarly, is uniquely defined by its purely outgoing 
behavior at spatial infinity 


€^( 0 . 




1/2 


P+2[p(a;)r* + [Mp^/p(a;)] ln(r/M)] 


(16a) 


and, at the horizon, it has an asymptotic behavior of the 
form 






(16b) 


In the previous expressions, the function p{uj) = 
{uP — denotes the “wave number” while the coef¬ 
ficients Bg~\uj) and B^^'^\uj) are com¬ 

plex amplitudes. Here, it is important to note that these 
coefficients as well as the in and up modes can be defined 
by analytic continuation on the two-sheet Riemann sur¬ 
face associated with p{uj). By evaluating the Wronskian 
W^(a;) at r* —>■ —oo and r* —> -boo, we obtain 

Wg{uj) = 2iu}Alr\u}) = 2iujB^/\co). (17) 
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D. Solution of the wave equation with source 


E. General expression of the source 


Using the Green’s function (13), we can show that the 
solution of the Regge-Wheeler equation with source (10) 
is given by 


f+°° / 2M \ 

(t>ujim{r) = J dr'^ Gc^£(r*,r') ( 1 - ) PuiimiK) 

(18) 

or, equivalently, by 

p+oo 

(t>u^im{r)= / dr'G^i{r,r')p^irn{r’)- (19) 

J 2 M 


As a consequence, the response ( 8 ) detected by an ob¬ 
server at (t, r, 9, ip) is given by 

-i-cxD +e 

<^{t,r,9,ip) = - EE 4>em.{t,r)Ye^(9,(p) (20) 

£—0 m——l 

where 

1 p+oo+ic 

(l)£m{t,r) = —— duje~''‘^*(t)ujem{r) ( 21 ) 

V 27r j _oo+ic 


(here c > 0 ) denotes the partial response corresponding 
to the (£, m) mode. 

For r —>■ -foo, the solution (19) reduces to the asymp¬ 
totic expression 


4^Lj£m (^ ) 


1 

\ui) 


UJ 

.pH. 


r+OO 


dr' pu,£m{r'). 


f2M 


( 22 ) 


This result is a consequence of Eqs. (13), (16a) and (17). 
However, it is very important to note that, in practice, 
if we consider that the observer is located at a finite dis¬ 
tance from the BH, it is not possible to use the asymp¬ 
totic behavior (16a) in Eq. (13) but it is necessary to 
obtain (^“^(r) numerically by solving carefully the Regge- 
Wheeler equation (14): this is the price to pay in order 
to take into account the dispersive nature of the massive 
scalar field. In that case, assuming that the observer is 
located beyond the source, we must replace ( 22 ) with 


4^Lj£'m{r) — 


1 


2iuiA^f, ^ {ui) 
/‘ + 00 


4>Zir) 


/2M 


dr'p^£^{r') (23) 


and then we have for the {i,m) waveform 


litA) = - 


V^J- 


^+ 00 + 2 C 


du! 


00+2C 
+ 00 


^ 2iwA^ ^ (w 


/ dr'(£^i{r')p^irn{r'). (24) 


/2M 


By extremization of the action (4), we can obtain the 
equations governing the timelike geodesic 7 followed by 
the massive particle. We shall denote by tp(r), ■rp(r), 
0 p(r) and Ppir) its coordinates and we remark that, 
without loss of generality, its trajectory can be consid¬ 
ered to lie in the BH equatorial plane, so we assume that 
0p(r) = tt/2. We have 


and 


/ 2M\ dtp _ E 

\ rp J dr mo’ 

2 dpp L 


(25a) 

(25b) 


dr J \ rp J rp 



2 

- 1 


(25c) 


where E and L are, respectively, the energy and the an¬ 
gular momentum of the particle. Here, it is important to 
recall that E and L are two conserved quantities. 

For the general trajectory previously described, the ex¬ 
pression ( 6 ) of the source is then given by 

p{x)=q5{9-TT/2) f -^Gdit-tpiT)) 

v /7 ‘p\‘) 

x6{r-rp{T))6{(p-(fpir)). (26) 


III. WAVEFORMS EMITTED BY A SCALAR 
POINT PARTICLE ON A PLUNGE 
TRAJECTORY 


A. Source due to a scalar point particle on a 
plunge trajectory 

In this section, we consider a particle plunging into the 
BH from the ISCO at nsco = 6 M or, more precisely, 
from slightly below the ISCO. In order to describe its 
trajectory, we need its angular momentum and its energy 
on the ISCO. They are determined by using Eqs. (Al) 
and (A2) which have been established for an arbitrary 
circular trajectory. We have 


-^ISCO 

mo 

= 2 v^A/ 

(27a) 

^^ISCO 

2^2 

(27b) 

mo 
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Substituting Lisco and Eisco into the geodesic equa¬ 
tions (25a)-(25c), we obtain after integration 


tpjr) 

2M 


2M 


-\-2 tanh 


(3M/r - 1/2)^/^ 


1 (6M/r - 1) 

2M 


1/2 

(28) 
























and 


0.03 


5 


, , 2-\/3r 

where to and ipo are two arbitrary integration constants. 
From (29), we can write the spatial trajectory of the 
plunging particle in the form 


This trajectory is displayed in Fig. 1. 


105* 76" 



FIG. 1. The plunge trajectory obtained from Eq. (30). Here, 
we assume that the particle starts at r = risco(l ~ e) with 
e = 10~^ and we take tpo = 0. The red dashed line at r = 6M 
and the red dot-dashed line at r = 2M represent the ISCO 
and the horizon, respectively. 

After integration and by using, in particular, the 
change of variable t — )• (fpir), Eq. (26) leads to 


p{t,r,e,^) = ^ 


r-hoo 


+CO +£ 


duj f 


m—— 




(6M - 


yl[\.q)\yuoA) ( 31 ) 


and we can then show that (31) can be written in the 
form (9) with 




3gyT gdwtp(r)-m(^p(r)] ^ , 


(32) 


B. Quadrupolar waveform produced by the 
plunging particle 

From Eq. (24) with the source term given by Eq. (32), 
we can now obtain numerically the partial waveform 
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FIG. 2. Quadrupolar waveforms produced by the plunging 
particle. The results are obtained for a massless scalar field 
{a = 0) and the observer is located at (a) r — lOM, (b) 
r = 20M and (c) r = SOM. 


4>emit,r) emitted by the plunging particle. In this sec¬ 
tion, we shall only focus on the (£ = 2, m = 2) mode of 
the scalar field and we shall emphasize the role of the 
mass parameter fi and the location of the observer. 

In Figs. 2, 3 and 4, we display the partial waveforms 
corresponding, respectively, to the values d = 0, 0.25 
and 0.35 of the reduced mass parameter. For each one, 
we consider that the observer is located at r = lOM, 
20M and 50M. The waveforms have been obtained by 
assuming that the particle starts at r = nsco(l ~ e) with 
e = 10“"^ and, furthermore, we have taken (/jq = 0 and 
to/(2M) = 500 in order to shift the interesting part of 
the signal in the window t/(2M) G [0, 600]. To construct 
these waveforms and, in particular, to obtain the func¬ 
tions and 4’'^^ as well as the coefficient we 

have numerically integrated the Regge-Wheeler equation 
(14) with the Runge-Kutta method. The initialization of 
the process has been achieved with Taylor series expan¬ 
sions converging near the horizon and we have compared 
the solutions to asymptotic expansions with ingoing and 
outgoing at spatial infinity that we have decoded by Fade 
summation. Moreover, in Eq. (24), we have discretized 
the integral over w. For d = 0, 0.25 and 0.35, in or¬ 
der to obtain numerically stable waveforms, we can limit 
the range of frequencies to —6 < 2Muj < -1-6 and take 
for the frequency resolution 2M5u) = 1/10000. It should 
be noted, however, that in order to extract correctly the 
spectral content of the signals (see Sec. HIE), it is nec¬ 
essary to work with higher-frequency resolutions which 
strongly depend on the mass parameter d. Eor example, 
for d = 0 it is sufficient to take 2M5u} = 1/10000, while 
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FIG. 3. Quadrupolar waveforms produced by the plung¬ 
ing particle. The results are obtained for a massive scalar 
field (3 = 0.25 is below the threshold 3^ ~ 0.2722) and the 
observer is located at (a) r = lOM, (b) r — 20M and (c) 
r = 50M. 
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FIG. 4. Quadrupolar waveforms produced by the plunging 
particle. The results are obtained for a massive scalar field 
(3 = 0.35 is above the threshold 3^ « 0.2722)and the observer 
is located at (a) r — lOM, (b) r = 20M and (c) r = SOM. 


for a = 0.35, we have worked with 2M6uj = I/IOOOOO. 

For the massless scalar field (see Fig. 2), the waveforms 
can be decomposed in three phases: (i) an “adiabatic 
phase” corresponding to the quasicircular motion of the 
particle near the ISCO (see Fig. 1), (ii) a ringdown phase 


due to the excitation of QNMs and (iii) a late-time phase. 
Such a decomposition remains roughly valid for the mas¬ 
sive field (compare Figs. 3 and 4 with Fig. 2) but now, as 
we shall see in more detail in the following subsections, 
the behavior of the signal is modified by the excitation of 
the QBSs and, moreover, at a large distance from the BH, 
it strongly depends on a threshold value of the dimen¬ 
sionless coupling constant a which separates the regions 
where the (£, to) mode of the field is propagating (but 
dispersive) or evanescent. This threshold corresponds to 
the mass parameter = 2Qigco where flisco denotes 
the angular velocity of the particle moving on the ISCO 
or, in other words, its orbital frequency. The existence 
of this threshold and its consequences are discussed in 
detail in Appendix A. 


C. The adiabatic phase and the circular motion of 
the particle on the ISCO 



FIG. 5. Quadrupolar waveforms produced by the plunging 
particle (blue line) and by a particle orbiting the BH on the 
ISGO (red dashed line). The reduced mass (3 = 0.25) is 
below the threshold (3c « 0.2722). 

In Figs. 5 and 6, we compare the quadrupolar wave¬ 
form produced by the plunging particle we have obtained 
in Sec. IIIB with the quadrupolar waveform produced 
by a particle orbiting the BH on the ISCO we discuss 
in Appendix A. We consider two particular masses given 
by a = 0.25 and a = 0.35, i.e. one below the thresh¬ 
old value die ~ 0.2722 and the other one above, but the 
discussion remains valid even for other masses if they 
are not too large. In the adiabatic phase, the waveform 
emitted by the plunging particle is described very accu¬ 
rately by the waveform emitted by the particle living on 
the ISCO. This can be easily understood by noting that 
the initial position of the plunging particle is very close 
to the ISCO so that it undergoes an adiabatic inspiral 
along a sequence of quasicircular orbits near the ISCO. 
As a consequence, the comments made in Appendix A 
(see also Fig. 13) permit us to understand the behavior 
of the waveform in the adiabatic phase and to interpret 
Figs. 3 and 4. It is, in particular, important to note that 
in the adiabatic phase: 
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— For a given “distance” r of the observer, the wave¬ 
form amplitude decreases as the mass increases and 
vanishes for large masses. 

— Above the threshold, for an observer at spatial in¬ 
finity, the waveform amplitude vanishes. (However, 
as we shall see in Sec. HIE, this result is slightly 
modified due to the excitation of the long-lived 
QBSs). 



t/(2M) 


FIG. 6. Quadrupolar waveforms produced by the plunging 
particle (blue line) and by a particle orbiting the BH on the 
ISCO (red dashed line). The reduced mass (S = 0.35) is 
above the threshold (oc ~ 0.2722). 


D. The ringdown phase and the excitation of 
QNMs 

In Figs. 7 and 8, we compare the quadrupolar wave¬ 
form produced by the plunging particle we have obtained 
in Sec. IIIB with the quadrupolar quasinormal waveform 
given by Eq. (B3) which corresponds to the fun¬ 
damental (€ = 2,n = 0) QNM we discuss in Appendix 
B. We have considered two locations for the observer 
(r = lOM and r = 50M) and three particular reduced 
masses given by ci = 0, 0.25 and 0.35. When the reduced 
mass a and the distance r are not too large, the quasi¬ 
normal waveform describes very accurately the ringdown 
phase. However, we have checked that if the reduced 
mass a or the distance r increases, the agreement is not 
so good (cf., e.g. Fig 8c where a and r are both large). 
In our opinion, this is due to the excitation of QBSs (see 
Sec. HI E) which blur the QNM contribution. 


E. Excitation of QBSs 
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FIG. 7. Gomparison of the quadrupolar waveform pro¬ 
duced by the plunging particle (blue line) and the quadrupo¬ 
lar quasinormal waveform (red dashed line). The results are 
obtained for an observer at r = lOM. 
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FIG. 8. Gomparison of the quadrupolar waveform pro¬ 
duced by the plunging particle (blue line) and the quadrupo¬ 
lar quasinormal waveform (red dashed line). The results are 
obtained for an observer at r = SOM. 


In this subsection, we focus on the spectral content of 
the waveform by considering separately the adiabatic and 
the late-time phases. The associated spectral contents 
can be obtained by using the Fourier transform machin¬ 
ery and limiting the time integrations to the considered 
phase. 


In Figs. 9 and 10, we display the spectral content of 
the adiabatic phase of the quadrupolar waveform pro¬ 
duced by the plunging particle. For a reduced mass pa¬ 
rameter d below the threshold dc (see Fig. 9), we can 
observe a peak at w = 2njsco- H is, of course, as¬ 
sociated with the quasicircular motion of the plunging 
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FIG. 9. Spectral content of the adiabatic phase of the 
quadrupolar waveform produced by the plunging particle. 
The results are obtained for a massive scalar field (S = 0.25 is 
below the threshold dc ~ 0.2722) and the observer is located 
at (a) r = 20M and (b) r = SOM. We observe the signature 
of the quasicircular motion of the plunging particle. 



2Maj 


FIG. 10. Spectral content of the adiabatic phase of the 
quadrupolar waveform produced by the plunging particle. 
The results are obtained for a massive scalar field {a = 0.35 is 
above the threshold dc « 0.2722) and the observer is located 
at (a) r = 20M and (b) r = SOM. We observe, in addition 
to the signature of the quasicircular motion of the plunging 
particle, that of the first long-lived QBS. 


particle near the IS CO. For a reduced mass parameter a 
above the threshold dtc (see Fig. 10), the same peak at 
w = 2D,isco is present. However, it is important to note 
that its height decreases very rapidly as the distance r 
of the observer increases. This is due to the evanescent 
character of the {£ = 2, m = 2) mode below etc- It is, 
moreover, very interesting to remark the presence of an¬ 
other peak at a frequency equals to the real part of the 
complex frequency of the first long-lived QBS (see Ta¬ 


ble II). In other terms, we can observe the excitation of 
the first QBS in the adiabatic phase! As we have already 
noted in Sec. IIIC, above the threshold, for an observer 
at a large distance from the BH, the contribution of the 
quadrupolar waveform associated with the quasicircular 
motion vanishes but the QBSs pick up the slack. 



2Mw 


FIG. 11. Spectral content of the late-time phase of the 
quadrupolar waveform produced by the plunging particle. 
The results are obtained for a massive scalar field (q = 0.25 is 
below the threshold ~ 0.2722) and the observer is located 
at (a) r = 20M and (b) r = 50M. We observe the signature 
of the first long-lived QBS. 



FIG. 12. Spectral content of the late-time phase of the 
quadrupolar waveform produced by the plunging particle. 
The results are obtained for a massive scalar field (q = 0.35 is 
above the threshold Qc « 0.2722) and the observer is located 
at (a) r = 20M and (b) r = 50M. We observe the signature 
of the first long-lived QBSs. 

In Figs. 11 and 12, we display the spectral content of 
the late-time phase of the quadrupolar waveform pro¬ 
duced by the plunging particle. We can observe the ex- 
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citation of QBSs (see also Table II). We can also note 
that, as the reduced mass parameter a increases, the 
spectrum of the frequencies of the QBSs spreads more 
and more and it is then possible to separate the different 
excitation frequencies. 


IV. CONCLUSION AND PERSPECTIVES 

If the graviton has a mass, the study of the gravita¬ 
tional radiation generated by a particle plunging into a 
BH is certainly a problem of fundamental importance 
whose solution could help us to test the various mas¬ 
sive gravity theories. Because this problem seems to us 
rather difficult, theoretically as well as numerically, we 
have chosen to work with a toy model where the massive 
spin-2 perturbations of the BH are replaced by a mas¬ 
sive scalar field and we have considered that this field 
and the plunging particle are linearly coupled. In our 
opinion, the study of the scalar radiation generated by 
the plunging particle has permitted us to exhibit some 
features that should also be present in the solution of 
the analogous problem in massive gravity because they 
are directly associated with characteristics shared by the 
massive scalar and spin-2 fields such as (i) the propagat¬ 
ing or evanescent nature of the partial modes and (ii) the 
existence of QBSs. 

In this article, by limiting our study to the case of the 
Schwarzschild BH and to the {£ = 2,m = 2) mode of the 
scalar field, we have more particularly shown that: 

— The waveform produced by a particle plunging into 
a Schwarzschild BH from slightly below the ISCO 
can be roughly decomposed in three phases: (i) an 
adiabatic phase corresponding to the quasicircular 
motion of the particle near the ISCO, (ii) a ring- 
down phase due to the excitation of QNMs and (iii) 
a late-time phase. 

— In the adiabatic phase, for an observer at a large 
distance from the BH, the behavior of the wave¬ 
form depends on a threshold value of the reduced 

mass parameter a. This threshold corresponds to 
the mass parameter /ic = 2Hisco where Hisco de¬ 
notes the orbital frequency of the particle moving 
on the ISCO. This threshold value separates the 
dispersive and the evanescent regimes. Above the 
threshold, for an observer at spatial infinity, the 
waveform amplitude vanishes. For a given distance 
r of the observer, the waveform amplitude decreases 
as the mass increases and vanishes for large masses. 

— The ringdown phase (oscillations and damping) is 
very well described from the excitation of the first 
QNM, i.e. the least damped one. 

— In the late-time phase, whatever the mass parame¬ 
ter, we can observe the excitation of QBSs. In the 
adiabatic phase, the excitation of QBSs only occurs 


for masses above the threshold. These QBSs could 
dominate the signal at a large distance from the 
BH. 

— For large masses, the ringdown phase is modified 
by the excitation of QBSs which blur the QNM 
contribution. 

It should be noted that, mutatis mutandis, the behav¬ 
iors of an arbitrary (£, m) waveform and of the quadrupo- 
lar waveform generated by the plunging particle are quite 
similar. In general, the waveform amplitude in the adia¬ 
batic phase decreases as £ increases or m decreases (see 
Figs. 14 and 15). However, it should be noted that this 
does not remain true around the threshold values for an 
observer at a large distance from the BH (see Fig. 14). 

We hope in the near future to extend this work to 
the massive spin-2 perturbations of the Schwarzschild BH 
and to consider higher values of the mass parameter in 
order to check if the intrinsic giant ringings predicted 
in our previous works [29, 30] are generated in physical 
processes. With this aim in view, there remain a lot of 
theoretical and numerical difficulties to overcome. 
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Appendix A: Waveforms produced by a scalar point 
particle living on the ISCO 

In this appendix, we provide a simple closed-form ex¬ 
pression for the emitted waveform when the particle lives 
on the ISCO and we analyze its behavior as the reduced 
mass parameter a increases and as the distance r of the 
observer changes. These results are helpful in order to 
describe the adiabatic phase of the waveform emitted by 
a point particle on a plunge trajectory (see Sec. HIC). 


1. Source due to a point particle on a circular orbit 
and associated waveform 

Here, we assume that the particle moves on a circular 
orbit with radius rg, i.e. we have rp(T) = rg = Const, 
in the geodesic equations (25). After integration, they 
provide for the angular momentum and the energy of the 
particle 


^ f Mrg \ 
mg \l-3M/rg) 


(Al) 
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and 


E _ (1 - 2M/ro) 

Wo (1 - SM/ro)^"'^' 


(A2) 


Moreover, we find that the angular coordinate tpp is given 

by 

if we use the proper time to describe the particle motion 
and by 


iPp{t) = nt (A-4) 

if we use the Schwarzschild time. In Eq. (A4), 


n = 



(A5) 


denotes the angular velocity of the particle (i.e. its orbital 
frequency). 

After integration and by using, in particular, the 
change of variable r —>■ ipp^r), we can show that (26) 
can be written in the form (9) with 


Puim ) 


ro \ ro J 

xS{r - ro)6{u} - mn)Yl^(J^,oY (A6) 


The (£, m) waveform produced by the particle orbiting 
the BH on a circular orbit with radius rg can now be ob¬ 
tained by substituting the source term (A6) in Eq. (24). 
After integration, we have 


4>tra{t,r) = 


q 

2iro 



1/2 

\Y*- 

^ im 





(A7) 


It should be noted that (i) 

can limit our study to the modes with m > 0, and (ii) 
(t’tra = 0 for ^ — TO odd. This will simplify the discussions 
in the next subsection. 


2. Waveforms for a particle on the ISCO 

We now specialize to the point particle living on the 
ISCO, i.e. we assume that rg = 6M in the waveform 
(A7) and we note also that (A5) reduces to 

We intend to highlight both the influence of the mass p of 
the scalar field and the distance r of the observer and to 



a 


FIG. 13. Waveform amplitudes for a particle on the ISCO 
and reduced masses in the range a £ [0,2]. The results are 
obtained from (A7) for the particular (f = 2, m = 2) mode. 
We study the role of the location of the observer. For an 
observer at large distance from the BH, the transition between 
the propagating (but dispersive) and evanescent behaviors of 
the response occurs at Sc ~ 0.2722. 


show more particularly that the classical and well-known 
results concerning the massless field cannot be naively 
extended to massive fields. 

In Fig. 13, we focus on the {I = 2,m = 2) wave¬ 
form. For an observer at spatial infinity, the behavior 
of the waveform amplitude depends in a complex way 
on the mass of the scalar field. This is due to the term 
p{uj = TOflisco) = [(’Tiflisco)^ ~ which is implic¬ 

itly present in the expression (A7) of the waveform. In¬ 
deed, depending on whether p < fj.^ ot p > Pc with 


/ic = toIIisco, (A9) 

this term is responsible for the propagating (but disper¬ 
sive) or evanescent character of the (£, to) mode of the 
massive field. In Fig. 13, we observe the decreasing of 
the waveform amplitude when the reduced mass param¬ 
eter a increases and we note that, above the threshold 
value etc corresponding to fic, the amplitude vanishes. 
As a consequence, above etc, the excitation of the system 
scalar field-BH by a particle orbiting on the ISCO can¬ 
not be observed at spatial infinity. This abrupt behavior 
remains valid for large distances (for r > SOM) but is 
smoothed for short distances with a vanishing of the am¬ 
plitude for large masses. As a consequence, for “very” 
massive scalar fields, the excitation of the system scalar 
field~BH by a particle orbiting ISCO cannot be observed 
whatever the distance. 

In Figs. 14 and 15, we show that the behavior of the 
(£ = 2, TO = 2) waveform can also be observed for ar¬ 
bitrary (i, to) waveforms. Moreover, in Fig. 14, it is 
interesting to note that, in a large domain around the 
particular threshold masses, the mode with the lowest £ 
does not provide the greatest waveform amplitudes. In 
this, a massive field is completely different from a mass¬ 
less one. Furthermore, in Fig. 15 we can observe that, for 
a given £, the waveform amplitude decreases with to and 
therefore that the {£, m = ±£) modes induce the high- 
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a 


FIG. 14. Waveform amplitudes for a particle on the ISCO 
and reduced masses in the range a € [0,2]. The results are ob¬ 
tained from (A7) for some {I, m = £) modes and an observer 
at (a) r = lOM, (b) r = SOM and (c) r = lOOM. The tran¬ 
sition between the propagating (but dispersive) and evanes¬ 
cent behaviors of the response occurs at otc ~ 0.2722 for the 
{£ = 2,m = 2) mode, at dc ~ 0.4082 for the (f = 3, m = 3) 
mode and at dc ~ 0.5443 for the (f = 4, m = 4) mode. 


est amplitudes. Such behavior is well known for massless 
fields and we can see that it remains valid for massive 
fields. 


Appendix B: Quasinormal frequencies, QNMs and 
associated ringings 

In this appendix, we consider the weakest damped 
QNM with £ = 2 and we obtain numerically its response 
when it is excited by the plunging particle. This result 
is helpful in order to describe the ringdown phase of the 
waveform emitted by a point particle on a plunge trajec¬ 
tory (see Sec. HID). 

The response of a QNM to an excitation (the so-called 
quasinormal waveform or, in other words the associated 
ringing) can be constructed from (i) the intrinsic charac¬ 
teristics of the BH (the quasinormal frequency toim the 
QNM itself and the corresponding excitation factor B^n 
- here n = 0 corresponds to the fundamental QNM, i.e. 
the least damped one and n = 1, 2 ,... to the overtones) 
and (ii) the source of the excitation. In the following 
subsections, we describe this construction. 
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a 

FIG. 15. Waveform amplitudes for a particle on the ISCO 
and reduced masses in the range a € [0,2]. The results are 
obtained from (A7) for modes with angular momentum £ = 4. 
We study the influence of the azimuthal number m for an 
observer at (a) r — lOM and (b) r = lOOM. 



1. Quasinormal frequencies and excitation factors 

We recall that the quasinormal frequencies uiin are the 
zeros of the Wronskian Wi{uj) given by Eq. (17) lying in 
the lower part of the first Riemann sheet associated with 
the function p(w) = (w^ — (see Fig. 16) and that 

the corresponding excitation factors are defined by 


Btn = 


1 

2p{uj) 


4+)(a;) 


dA 


i _V 

du 


(Bl) 


It should be noted that when the Wronskian Wi{iS) van¬ 
ishes (i.e. if w = Loin), the functions and (j)^^ are 
linearly dependent and propagate inward at the horizon 
and outward at spatial infinity, such behavior defining 
the QNMs. 

We also recall that the complex spectrum of the quasi¬ 
normal frequencies is symmetric with respect to the imag¬ 
inary w axis. In other words, if is a quasinormal 
frequency lying in the fourth quadrant, is the sym¬ 

metric quasinormal frequency lying in the third one. In 
Table I, we have considered the least damped £ = 2 QNM 
and we have given its quasinormal frequency a;2o with 
positive real part and the associated excitation factor ,620 
for three particular values of the reduced mass parameter 
a. The results have been obtained by using the numerical 
methods described in Sec. II of Ref. [30]. 


2. Quasinormal waveforms 

We now deform the contour of integration in Eq. (24) 
in order to extract a residue series over the quasinormal 
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TABLE I. Evolution of the quasinormal frequency and the 
excitation factor of the fundamental 1 = 2 QNM. 


(£, n) a 2 MuJin Bin 

(2,0) 0 0.96729 - 0.193521 0.11935 + 0.01343i 

(2,0) 0.25 0.97717 - 0.19012i 0.12241 + 0.00737i 


(2,0) 0.35 0.98669 - 0.186841 0.12548 + 0.00091i 





FIG. 16. The two Riemann sheets of the function p{uj) = 
{uP — and the branch cut we have chosen. We have in¬ 

troduced the polar coordinates {pi,9i) and (^ 2 ,^ 2 ) such that 
(jj — fj, = and ui + p = p 2 e*®^. Here 0 < 61,62 < 2%. 


contribution of i.e. its symmetric with respect to 

the imaginary axis. The expression (B3) has been sim¬ 
plified by using some symmetry properties in the change 
f particular, 


P{-^L) = -[Pi^^in)]*, (B6a) 

(B6b) 

= iZU*- (B6c) 


It should be noted that, in our problem, the spherical 
symmetry of the Schwarzschild BH is broken due to the 
asymmetric plunging trajectory. It is this dissymmetry 
which, in connection with the presence of the azimuthal 
number m, forbids us to gather the two terms in Eq. (B3). 

Let us hnally remark that does not pro¬ 

vide physically relevant results at “early times” due to 
its exponentially divergent behavior as t decreases. It 
is necessary to determine, from physical considerations 
when this is possible or by using a numerical approach, 
the time beyond which this waveform can be used, i.e. 
the starting time <start of the BH ringing. 


Appendix C: Complex frequencies of the first QBSs 


frequencies. It is given by 


+00 

= (B2) 

n—0 

with 

=V2Tr (Cemn 

+ 25^-n[C„^(r*)]*e--*). (B3) 

Here Cimn and Vimn denote the extrinsic excitation co¬ 
efficients defined by 


Ci 


mn 


J 2 M [a;/p(a;)] A^^\ijj) 


and 


(B4) 




I 2 M [uj/p{uj)] A'y^\u}) 


(B5) 


The first term in Eq. (B3) is the contribution of the quasi¬ 
normal frequency lying in the fourth quadrant of the 
first Riemann sheet of p(a;) while the second one is the 


The complex frequencies of the QBSs are the zeros 
of the Wronskian Wi{uj) given by Eq. (17) lying in the 
lower part of the second Riemann sheet associated with 
the functionp(w) = (w^—(see Fig. 16). Their spec¬ 
trum is symmetric with respect to the imaginary w axis. 
They can be numerically obtained by using the method 
which has permitted us to determine in Sec. B 1 the quasi¬ 
normal frequencies but now, because we are working on 
the second Riemann sheet associated with p{ui), it is nec¬ 
essary to use —p{oj) instead of p{uj) (see also Fig. 16). 
In Table H, we have given a sample of the complex fre¬ 
quencies of the long-lived QBSs relevant to the spectral 
content of the waveform emitted by the plunging particle 
(see Sec. HIE). 


TABLE 11. Complex frequencies uJe„ of the first long-lived 
QBSs. 


(f,n) 

a 



(2,0) 

0 


— 

(2,0) 

(2,1) 

0.25 

0.24978 

0.24988 

- 6.08009 X 10-^'^i 

- 3.63454 X 10-^'^i 


(2,0) 0.35 0.34940 - 1.08832 x 10" 

(2.1) 0.34966 - 1.03797 x 10" 

(2.2) 0.34978 - 3.90581 x 10' 

(2.3) 0.34985 - 2.42608 x 10' 

(2.4) 0.34989 - 1.59115 x 10' 
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